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Abstract 

We perform model independent and model dependent analyses of solar neutrino 
data including the neutral current event rate from SNO. The inclusion of the first SNO 
NC data in the model independent analysis determines the allowed ranges of S B flux 
normalisation and the v e survival probability more precisely than what was possible 
from the SK and SNO CC combination. We perform global z/ e — inactive oscillation 
analyses of solar neutrino data using the NC rate instead of the SSM prediction for 
the 8 B flux, in view of the large uncertainty in the latter. The LMA gives the best 
solution, while the LOW solution is allowed only at the 3a level. 



1 e-mail: abhi@theory.saha.ernet.in 
2 email : sandhya® hep . phys .soton.ac.uk 
3 e-mail: sruba@mri.ernet.in 
4 e-mail : dproy @ theory, tifr . res .in 

1 



The neutral current results from the Sudbury Neutrino Observatory measures for the 
first time the total flux of 8 B neutrinos coming from the Sun [Q . In a recent paper || we had 
examined the role of the anticipated NC data from SNO in enhancing our understanding 
of the solar neutrino problem. The SNO NC rate can be expressed in terms of SNO CC 
and SK elastic scattering rates as || 

r>NC _ r>CC , / r>el r>CC w /,\ 

SNO SNO ' \ SK n-SNOJl'i \ L ) 

where r = a^ T /a^ +NC ~ 0.157 for a threshold energy of 5 MeV (including the detector 
resolutions and the radiative corrections to v — e scattering cross-sections). All the rates 
are defined with respect to the BBP2000 Standard Solar Model (SSM) [|J. We showed in 
@ that because SNO has a greater sensitivity to the NC scattering rate as compared to 
SK, the SNO NC measurement will be more precise and hence incorporation of this can be 
more predictive than the SNO CC and SK combination. We took three representative NC 
rates - R S N N C ° = 0.8,1.0 and 1.2 (±0.08) and showed that 

1. For a general transition of v e into a mixture of active and sterile neutrinos the size of 
the sterile component can be better constrained than before. 

2. For transition to a purely active neutrino the 8 B neutrino flux normalisation and the 
survival probability P ee are determined more precisely. 

3. We had also performed global two flavour oscillation analysis of the solar neutrino 
data for the v e — v active case, where instead of Rsk and R$no we use d the quantities 
Rsk/Rsno an d R-sno/ -R-sno- These ratios are independent of the 8 B flux normalisa- 
tion and hence of the SSM uncertainty. We showed that use of these ratios can result 
in drastic reduction of the allowed parameter regions specially in the LOW-QVO area 
depending on the value of the NC rate. 

We now have the actual experimental result 

Rggo = 1.01 ± 0.12 (2) 

while eq. ([I]) gives 1.05 ±0.15. Thus in 306 live days (577 days) the SNO NC measurement 
has achieved a precision, which is already better than that obtained from the SK and SNO 
CC combination. This paper follows closely the analysis that we have done in [Q] but 
incorporating the actual data. In addition we also perform an alternative global analysis 
for v e — v ac ti ve oscillation by letting the 8 B normalisation factor fs vary freely, where the 
inclusion of Rsno( = Sb) m the At serves to control this parameter. As we shall see below 
the two methods of global analysis give very similar results. 

In section 1 we discuss the constraints on the electron neutrino survival probability, 
the 8 B normalisation factor fs and the fraction of sterile component without assuming any 
particular model for the probabilities. In section 2 we perform the global analyses assuming 
two flavour v e — v act i ve oscillation. 
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1 Model Independent Analysis 



For the general case of v e transition into a combination of v a ctive (^a) and v ste riie {y&) states 
one can write the SK, SNO CC and SNO NC rates as 

f B Pee + f B rP ea , (3) 
f B Pee, (4) 
f B (Pee + Pea), (5) 

where P ee and P ea denote the probabilities folded with the detector response function and 
averaged over energy. To extract a model independent bound on P ee one has to ensure an 
equality of the response functions which amounts to slight adjustment of the SK threshold 
energy and the rate || ||. Our approach is slightly different. We treat P ee to be effectively 
energy independent. The SK spectrum data indicates a flat probability down to 5 MeV 
0. This is corroborated by SNO |8|, |9| which now has a threshold of 5 MeV for kinetic 
energy of the observed electron. Hence we consider this assumption as justified and expect 
the results to be insensitive to the differences in the response functions. It should be noted 
here that in cotrast to the SNO CC events their NC events correspond to a neutrino energy 
threshold of 2.2 MeV. However it is clear from eq. (5) that for a u e to u a transition there 
is no reason to expect any energy dependence in Rg^ Q . On the other hand for the general 
case of v e transition into a combination of v a and v s our approach effectively assumes P es 
to be energy independent down to 2.2 MeV. A comparison of the current values R^no w hh 
P-sno * s shown in fig. 1. It constitutes a 5.3a signal for transition to a state containing an 
active neutrino component or alternatively a 5.3a signal against a pure sterile solution. 

Next we consider the general case where v e goes to a mixed state = v a sin a + u s cos a. 
Then one can write P ea = sin 2 a(l — P ee ). Substituting this in the equations (||) and (|5|) 
and eliminating P ee using equation (£|) one gets the following set of equations for fs and 
sin 2 a [| 

sinM/B-i^o) = {Rf K -R C s C Nc)/r, (6) 
sin 2 a(f B - R C sno) = r sno ~ r sno- ( 7 ) 

We treat sin 2 a as a model parameter. And for different input values of sin 2 a we determine 
the central value and the la and 2a ranges of fs by taking a weighted average of the 
equations (||) and (^). The corresponding curves are presented in fig. 2. Combining the 2 
a lower limit of fs from this fit with the 2a upper limit from the SSM (vertical lines) gives 
a lower limit of sin 2 a > 0.45 i.e. the probability of the active component is > 45%. Note 
that there is no upper limit on this quantity - i.e. the data is perfectly compatible with v e 
transition into purely active neutrinos. 

Assuming transition into purely active neutrinos (P ea = 1 — P ee ) we show in fig. 3 
the la and 2a contours in the Jb — Pee plane from the combinations SK+SNOCC and 
SK+SNOCC+SNONC. The inclusion of the NC rate narrows down the ranges of Jb and 
P ee . The error in after the inclusion of NC data is about half the size of the corresponding 
error from SSM as is seen from fig. 3. 
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2 Model dependent analysis 



In this section we present the results of our x 2 analysis of solar neutrino rates and SK 
spectrum data in the framework of two flavour oscillation of v e to an active flavour. We use 



the standard techniques described in our earlier papers |10|, JTJ excepting for the fact that 
instead of the quantities Rij K and Rsno we now 

fit the ratios Rsk/ r sno and r sno/ r sno- 
The 8 B flux normalisation gets cancelled from these ratios and the analysis becomes in- 
dependent of the large (16-20%) SSM uncertainty associated with this. We include in our 
global analyses the 1496 day SK zenith angle spectra jl5| . Since we use both SK rate and 
SK spectrum data we keep a free normalisation factor for the SK spectrum. This amounts 
to taking the information on total rates from the SK rates data and the information of 
the spectral shape from the SK spectrum data. The SNO CC and NC rates have a large 
anticorrelation. We have taken into account this correlation between the measured SNO 
rates in our global analyses. Further details of this fitting method can be found in @. In 
Table 2 we present the best-fit parameters, Xmin an d goodness of fit (GOF). The best-fit 



comes in the HIGH(LMA) region as before [11, 12]. However as is seen from fig. 4a the 



incorporation of the NC data narrows down the allowed regions, and in particular the LOW 
region becomes much smaller. 

We have also performed an alternative \ 2 fit to the rates of Table 1 ]l], |7|? HH along 



with the 1496 day SK spectra keeping f B as a free parameter. Even though we allow 
fs to vary freely the NC data serves to control within a range determined by its error. 
As we see from Table 2 and fig. 4b the results of this fit are very similar to the previous 
case. The best fit comes from the HIGH(LMA) region, while no allowed region is obtained 
for the LOW solution at the 99% CL level. Maximal mixing is seen to be disallowed at 
the 3a level. To illustrate the impact of the NC rate on the oscillation solutions we have 
repeated the free fg fit without this rate. The results are shown in fig 4c. Evidently the 
NC data plays a pivotal role in constraining the oscillation solutions, particularly in the 
LOW/QVO region, which is allowed only at the 3a level. It puts an upper bound on the 
Am 2 in the LMA region and rules out maximal mixing. 



3 Summary and Conclusions 

The first SNO NC data constitutes a 5.3a signal for transition into a state containing an 
active neutrino component. The inclusion of this data puts much tighter constraints on 
fs and P ee from a model independent analysis involving active neutrinos as compared to 
the SNO CC/SK combination. In this paper we have discussed two useful strategies, of 
incorporating the NC data in the global \ 2 analysis of rates and spectrum data, by which 
one can avoid the large 8 B flux uncertainty from the SSM. 

• We fit the ratios of the SK elastic and SNO CC rates w.r.t the NC rate, from which the 
fs cancels out. 

• We fit the rates by keeping fs as a free parameter, where the inclusion of the SNO NC 
rate (= fs) serves to control this parameter. 

Both the analyses give very similar results. They clearly favour the HIGH (LMA) so- 
lution, while a limited region of the LOW solution is also acceptable at the 3a level. The 
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maximal mixing solution is disfavoured at the 3a level. As more data accumulate one 
expects a substantial reduction in the error bar of the SNO NC rate, resulting in further 
tightening of the allowed regions of neutrino mass and mixing. 



Note Added: The paper [16] appeared on the net after completion of our work. In the 



region of overlap our results agree with theirs as well as with the updated version of |17] . 
It may be added here that the SNO CC and NC rates given in Table 1 are obtained as- 
suming undistorted energy spectra above 5 MeV, which for transitions to active neutrinos 
has good empirical justification as mentioned above. We thank Prof. Mark Chen of SNO 
collaboration for communication on this point. 
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experiment 


R 


composition 


Ga 


0.553 ± 0.034 


pp(55%),Be(25%),B(10%) 


CI 


0.337 ± 0.030 


B(75%),Be(15%) 


SK 


0.465 ± 0.014 (0.363 ± 0.014) 


£(100%) 


SNO{CC) 


0.349 ± 0.021 


5(100%) 


SNO{NC) 


1.008 ± 0.123 


5(100%) 



Table 1: The observed solar neutrino rates relative to the SSM predictions (BP2000) are 
shown along with their compositions for different experiments. For the SK experiment the 
v e contribution to the rate R is shown in parantheses assuming v e — > v a transition. In the 
combined Ga rate we have included the latest data from SAGE and GNO. 



Data 


Nature of 


Am 2 




tan 2 # 


Xmin 


Goodness 


Used 


Solution 


in eV 2 








of fit 


Ga + 


LMA 


9.66 x 10" 


-5 


0.41 


35.95 


80.08% 


SK/NC + 


LOW 


1.04 x 10" 


-7 


0.61 


46.73 


36.09% 


CC/NC + 


VO 


4.48 x 10" 


10 


0.99 


54.25 


13.84% 


SKspec 


SMA 


6.66 x 10" 


-6 


1.35 x 10~ 3 


67.06 


1.41% 


CI + Ga + 


LMA 


6.07 x 10" 


-5 


0.41 


40.57 


65.99% 


SK + CC + 


LOW 


1.02 x 10" 


-7 


0.60 


50.62 


26.14% 


NC + SKspec 


VO 


4.43 x 10- 


10 


1.1 


56.11 


12.39% 


+ Jb free 


SMA 


5.05 x 10" 


-6 


1.68 x 10" 3 


70.97 


0.81% 



Table 2: The Xmim the goodness of fit and the best-fit values of the oscillation parameters 
obtained for the analysis of the global solar neutrino data. 
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Figure 1: The SNO CC and NC rates shown relative to their SSM predictions. The 
dashed line is the prediction of the pure v e to v s transition. The pure sterile solution is 
seen to be disfavored at 5.3o\ 



8 




Figure 2: Best fit value of the 8 B neutrino flux fg shown along with its la and 2a limits 
against the model parameter sin 2 a, representing v e transition into a mixed state [y a sin a + 
z/gCosa). The dashed line denote the ±2<r limits of the SSM. 
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Figure 3: The la and 2a contours of solutions to the 8 B neutrino flux and the v e survival 
probability P ee assuming v e to v a transition. The la SSM error bar for Jb is indicated on 
the right. 
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Figure 4: The u e — > v a oscillation solutions to the global solar neutrino data using (a) Ga 
rate, the SK zenith angle energy spectra and the SK and SNO(CC) rates, both normalised 
to the SNO(NC) rate and (b) total Ga, CI, SK, SNO(CC) and SNO(NC) rates along 
with the SK zenith angle energy spectra, keeping the 8 B flux normalisation free. In 
both cases we use the SNO(NC) error as the error in the 8 B flux. The case(c) is similar 
to (b), but without using the SNO(NC) rate. 
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